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hadrons. The efficiency of the operators is shown by looking at the wave function of the first excited state, which 
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1. Introduction 

In quenched lattice gauge theory, hadronic 
states are created by acting with light quark (or 
anti-quark) creation operators on the vacuum. 
The ground state mass is then determined in the 
large (euclidean) time limit. In practice, however, 
only a limited range of t is available, so methods 
are needed to create the hadronic ground state 
efficiently. 

Since the main contamination of the ground 
state signal at larger t values comes from the first 
excited state, one requires a hadronic creation op- 
erator which maximises the ground state relative 
to this first excited state. 

One clear motivation for a trial hadronic opera- 
tor comes from considering heavy quark hadrons: 
for heavy mesons {cc and bb) the adiabatic ap- 
proximation is well justified, so that they can 
be modelled as heavy point particles (quarks) 
bound by a central potential V{R) between static 
colour sources. On a lattice efficient gluonic op- 
erators which create such a colour flux between 
static sources can be constructed using an itera- 
tive fuzzing algorithm |l|. This prescription cre- 
ates gluonic fluxes with a dominant ground state 
and very few excited states. To extend this ap- 
proach to lighter quarks, we use a lattice con- 
struction of a colour flux tube of length R to join 
two light quarks in a gauge invariant manner. By 
varying R we can explore the relative amplitude 
of ground state and excited state hadron created. 

This relative amplitude is usually called the 
Bethe Saltpeter (BS) wave function of the hadron. 
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Figure 1. The effective mass for the p (in lattice 
units) using the effective operators LL (-I-), and 
LF with R=4 (x), 8 (o), 10 (*) and 12 (D). 



It is the overlap between a quark and antiquark 
at distance R apart and the hadronic state which 
is an eigenstate of the hamiltonian (transfer ma- 
trix on a lattice). We use a fuzzed gluon flux 
prescription introduced in an earlier work with 
a similar approach y] to join the quarks, which 
corresponds to the adiabatic wavefunction as de- 
fined by Q. Some previous work has used quark 
and antiquarks in the (spatial) coulomb gauge in- 
stead [g[|| . This is less efficient (in our sense) and 
also leads to problems with image sources in the 
spatial periodic boundary conditions. 
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Figure 2. The effective mass for the Nucleon. 
Here we show the LLL, LLF and LFF operators 
for the various R values (see text). 



2. Lattice Measurements 

We use hght quark propagators in the 24^ x 48 
configurations at (3 = 6.2 and K — 0.f4144 us- 
ing the clover action obtained by UKQCD yj. 
Our most comprehensive wave function results 
come from an analysis of 12 configurations, al- 
though some quantities are available from a larger 
sample of 60 configurations. At these values of 
(3 and K the inverse lattice spacing is a^^ — 
2.73(5)GeV (determined from the string tension), 
while rriT^/vip « 0.77. 

We construct fuzzed gluon flux tubes following 
the successful methods used for studying the po- 
tential 0]: 
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Varying the parameters of the fuzzing prescrip- 
tion in our case gives relatively little change. 
Hence we select a smaller fuzzing level (5) and 
coarser fuzzing (c — 2) to minimise computer re- 
sources. 

Since we have access only to light quark propa- 
gators from a single source (0,0), we use a conven- 
tional local hadronic operator at the source. At 



the sink at time i, we use the spatially extended 
operator of length R along a lattice axis. A sum 
over spatial positions (to have momentum zero) 
and a sum over all 6 orientations (to get correct 
J^'-') is used. 

Here we follow the definitions and construction 
of the hadronic interpolating fields used by e.g. 
UKQCD [Q. The resulting hadronic correlation 
averaged over all configurations is then fitted to 
a two exponential function 

< h{0)H^{t, i?) > = co(-R)e-™°' + ci(i?)e-™i* 

We use a simultaneous fit to data at all R values 
and t values by making use of modified correlated 
X^ fits, which models the correlations between dif- 
ferent data points ^ . 

In the following we use R — 4, 8, 10 and 12, 
while the local operator corresponds to i? = 0. 
It will be shown below that this range of R val- 
ues is broad enough for the hadrons considered 
here. Statistical errors are determined using a 
bootstrap analysis. 

For the mesons we form local-fuzzed (LF) cor- 
relations by replacing one local propagator in the 
usual (LL) formalism by the fuzzed propagator 
discussed above. The effective mass results for 
vector meson are shown in Fig. 1. There are five 
observables corresponding to LL, and LF with the 
four different R values stated above. 

For the nucleon, we consider two different non- 
local operators. These involve a di-quark sepa- 
rated from a quark (single fuzzed - LLF) and an 
arrangement of three quarks all separated (dou- 
ble fuzzed - LFF). Fig. 2 shows the results thus 
obtained for the nucleon. 

Both figures clearly show that, by using fuzzed 
non-local operators, the plateau in the data starts 
at smaller values of t than for the purely local 
observable. The fuzzed data have larger errors at 
lower t (compared to the unfuzzed ones), but for 
larger t the errors are comparable. The unfuzzed 
contribution thus has an accurate but irrelevant 
component. 

The two exponential fit, which is stabilised by 
using several hadronic correlations to fit simulta- 
neously, has the advantage that the ground state 
can be better isolated than by using a single ex- 
ponential fit. 
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Figure 3. The tt wave function for the ground 
state (x) and first excited state (o), as functions 
of the length of the gluonic flux tube. 



Figure 4. As in Fig. 3, but for the Nucleon, using 
the LLF effective operators. 



3. Wave Functions 

To obtain a reasonably accurate determination 
of the wave function of the first excited state 
with our limited statistics, we fix the difference 
between the ground state and first excited state 
masses for the hadrons from results obtained re- 
cently by the UKQCD Collaboration ||l. These 
consist of hadronic propagators smeared at the 
source and sink (SS) and at the source only 
(SL) by applying the Jacobi smearing method on 
the existing 60 configurations at the origin, ob- 
tained on the 24^ x 48 lattice at /3 = 6.2 with 
K = 0.14144 [Q. A big advantage of having SS, 
SL and LL operators is that it allows a factoris- 
ing fit. These in turn provide tight constraints on 
the ground and excited state masses . We then 
proceed as before, fitting the local and non-local 
quark propagators calculated from the subset of 
12 configurations to the fit function given above. 

The wave functions for the ground and first ex- 
cited states obtained from the local operators and 
those involving only one fuzzed link, normalised 
at distance zero, are shown in Figs. 3 and 4 for 
the TT and Nucleon respectively. The wave func- 
tions for the p and A show similar behaviour. The 
results for the ground state wave function are in 
agreement with those obtained previously in the 
literature using similar gauge invariant definitions 



of the wave function 

As far as we are aware, the excited state wave 
function has not been studied either for mesons or 
baryons. The interesting feature that can clearly 
be seen for all the hadronic observables consid- 
ered here is the presence of a node in the excited 
state wave function as a function of the length 
of fuzzed links connecting the propagator ends at 
the sink, where the ratio of ground state wf to 
excited state wf becomes zero. This particular 
spatial extent (i? « 8 or « 3GeV~^ in physical 
units), which is more or less the same for all the 
hadronic observables, thus seems to be an optimal 
choice for producing a clean ground state mass, 
since the contamination of the ground state by 
higher excited states has been minimised. 
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